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Based on local unitary operators acting on a n- dimensional 
Hilbert-space, we investigate selective and collective operator 
basis sets for TV-particle quantum networks. Selective cluster 
operators are used to derive the properties of general cat- 
states for any n and N . Collective operators are conveniently 
used to account for permutation symmetry: The respective 
Hilbert-space dimension is then only polynomial in TV and 
governed by strong selection rules. These selection rules can 
be exploited for the design of decoherence-free subspaces as 
well as for the implementation of efficient routes to entangle- 
ment if suspended switching between states of different sym- 
metry classes could be realized. 



I. INTRODUCTION 

Problem-adapted representations are convenient tools 
for dealing with concrete physical models in virtually any 
branch of physics: Examples are the choice of coordinates 
in classical mechanics, of mode representations in linear 
wave theory, of state vector- or matrix-representations 
and of complete operator basis sets in quantum mechan- 
ics. The mode of adaption may refer to the internal sym- 
metry of the system under consideration and/or its cou- 
pling to the outside world (means of measurement and 
control). Adapted representations - though in principle 
equivalent to any other - are expected to simplify numer- 
ical calculations and to enhance insight, (cf. Q) 

Operators in quantum mechanics may represent ob- 
servables, states or transformations. In any case it is 
convenient to think in terms of "elements", i.e. basic op- 
erator sets, out of which any other operator could be con- 
structed : For a n-dimensional Hilbert-space there are 
n 2 such (orthogonal) basis operators (defining Liouville- 
space) . If we prefer to think in terms of basic observables 
or states, these basis operators should be chosen hermi- 
tian; a pertinent example are the 5'[/(n)-generators. If 
we rather think in terms of basic " actions" , the basis- 
operators should be unitary (thus defining basic unitary 
transformation). The latter approach has become the 
method of choice for investigations relating to quantum 
computation and quantum information processing [9: 
there we are typically concerned with sequences of uni- 
tary transformations. (Only for n — 2 are unitarity and 
hermiticity compatible requirements for a complete basis 
set.) 



Quantum networks (composite systems) may be de- 
scribed in terms of product-operators, Q; if each local op- 
erator is taken from the respective unitary basis set, also 
the product-operators are unitary. Furthermore, they are 
completely specified by the type of " action" to be applied 
on each subsystem n, fi — 1, 2, . . . , N; these Q-operators 
will be termed " selective" . Alternatively, we may intro- 
duce operators E which specify the action but not the 
" address" . In this case we are naturally led to " collective 
operators", defining a specific action on a given number 
a of subsystems. Complete sets require the inclusion of 
phases. Collective as well as selective operator sets are 
equivalent: in particular, we can express one type by 
the other. A subset of collective operators has permuta- 
tion symmetry: These are the only allowed operators for 
fundamentally indistinguishable subystems (fermions or 
bosons). 

Typical scenarios realized, e.g. in nanostructures, 
will exhibit neither complete selectivity nor complete 
non-selectivity (permutation symmetry). However, sys- 
tems with weakly restricted selectivity should still more 
efficiently be described by Q-Operators, systems with 
weakly broken permutation symmetry by collective op- 
erators. It is this latter theme of operational (partial) 
indistinguishability which will be of central interest in 
our present investigation. 

Our paper is organized as follows: In section [TJ we 
discuss local basis operators with special emphasis on 
unitary operators and introduce the complex coherence 
vector as a generalization to the Bloch vector. Section HI 
extends this concept to quantum networks by the use of 
index-selective cluster operators. Alternatively, section 
[V proposes the use of collective operators for which all 
subsystems are treated on equal footing. Applications of 
these concepts are worked out in section |y|, starting with 
highly entangled states and generalized cat states. Per- 
mutation symmetry plays the key role in the remaining 
applications showing the fundamental difference between 
selective and collective treatment of quantum networks. 



II. BASIS OPERATORS 

A. Transition- and SU(n)-operators 

The Hilbert-space of dimension n is taken to be 
spanned by a complete and orthonormal set of states 
j = 0, 1, 2, . . . , n— 1. They may be considered eigenstates 



1 



of some operator A, 



A 



3=0 



(2.1) 



where P\j — \i) (j \ ; = Pji denote n 2 transition op- 

erators, which are characterized by two state indices (= 
quantum numbers = eigenvalues Aj); they are orthonor- 

malized, tr |-Py-fjy} = ^ik^jh (tr = trace operation) and 

complete. Any other operator B s can then be written as 



B, 



7 j BsAjPij 



(2.2) 



where B Sji j = tr|s s ^|. We require the B s: (s 

.2 



0, 1, . . . , n z — 1) to form a complete orthogonal set, nor- 
malized to n, 



tr{s s 4t,| = ^ Ss .. B * 



(2.3) 



(Remember that the identity operator 1 is also normal- 
ized to 7i.) Note, that the transformation eq. (2.2) 
does not change the representation of any given op- 
erator; rather we change the set of operators within 
the ^-representation (i.e. the set of elementary ma- 
trices). A well-known example are the n 2 hermitian 



SU(n)-gencrators A = {A s 

{1,U X,Uq2, ■ ■ ■ ,U12, ■ ■ ■ ,#01, 



= 0,1,..., n 2 
,w ,.. .,W n - 2 } 



Uik 



Wi = 



n-l 
3=0 



■Phi) 



-i{P ik -Pu) 



-(Po> 



1} = 

(2.4) 

(2.5) 
(2.6) 
(2.7) 



(f+l)(Z + 2) 
+P H - (l + l)P l+u+1 ) , 

with tr|A s A S '| = n5 StS > for all s, s'. (The SU(n)- 

generators are often normalized to 2 rather than n 
the latter choice is more convenient for our purposes.) 
For n = 2 the representation of the Aj = Ao,Woi,Uoi, wq 
leads to the well-known Pauli matrices. For n > 2 the 
SU(n)-algebra tends to loose its convenience, because the 
corresponding operators are no longer unitary and their 
definition becomes rather unwieldy (cf. 0). 

With the unitary operators to be discussed next, the 
SU(n)-operators share two properties: 



i. The set consists of 



1 traceless operators B s , 



s^O and the unit-operator B = 1. 



ii. The operators based on projection operators Pjj 
are kept separated from the others. 

As a result of L . the p roduct-operators of a composite 
system (cf. Sect. p|jV| can be decomposed into a hier- 
achy of m-cluster-operators Q, where m is the number of 
indices s^O. 

As a result of ii., and if A = H, the Hamiltonian of 
the system under consideration, the expectation values of 
the Pjj or combinations of those are constants of motion 
(under the unitary evolution generated by H). 



B. Unitary operators 

The operators §,§ 
n-l 

U ab := V^ fc | fc + q) (k\ ; uj n ■- 

k=0 



(2. 



are defined in a double index notation a,b € [0, n— 1], 
where addition modulo n is denoted by underlining and 
the constant u> n is the n-th root of unity. We note for 
later reference, that 



]T ^ = nS ob ; E w » ^ = nSkk ° > 



(2.9) 



k=0 



b=0 



where 5y is the Kronecker-delta. The use of two in- 
dices leads to an easy interpretation of the action such 
an operator has on a state \k), 



Uab \k) = < fe Ifc + o) 



(2.10) 



namely b induces a phase shift and a causes a state-shift 
(for a concrete example see Appendix A). This set of n 2 
operators defines an orthonormal and complete operator 
basis for the Liouville-space with 



tr 



{UabUl} 



nS ac S t 



A operator A can be expanded like 

n-l 



1 

A = - Y] u ab U ab ; 



u ab := tr 



a,fc=0 



{ulA} 



(2.11) 



(2.12) 



Coming next we will state some of the basic properties 
these unitary operators obey. Because Uqq = 1, the or- 

thonormality relation ( |2.11 ) implies tr |?7 ab | = for 

a, b 0. This tracelessness will be essential for the def- 
inition of cluster-operators (cf. Sect. III). The adjoint 
operators are again, barring a phase, members of the set, 



U. 



u: 



■ U) n U- a -b ■ 



(2.13) 
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A very useful property of these operators is their cyclic 
symmetry, implying that any product reduces to just one 
operator of the set: 



UabUcd = ^nU a +c,b+d 



U ab U, 



, ,cd , — bcTj 
LO n U n Ua-cb-d 



U cd U a b = LO^UJ n a Ua-c.b-d 



UabUcdUef = u; b n c J n b+d ^U a 



±c+e.b+d+f 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



It is remarkable to note that the cyclic property allows 
to use - instead of the complete set {U a b} - the two op- 
erators {Z7o,n-l> ^71-1,0} only- All the others can then 
be generated as specific product forms, e.g. for n = 3: 

^00 = (U20) ,Uw = (U20) ,^22 = U20 ■ Uq2,Uqi — 

{Uo2) 2 ,Ui2 = (U 2 q) 2 ■ t/02, etc. This property is remi- 
nescent of the creation- and destruction operators, aJ , a, 
conveniently applied to harmonic oscillator models; it in- 
troduces a kind of " non-linearity" if expectation values of 
products (cf. eq. (2.28)) are approximated by products 
of expectation values. 

From eq. ( fTj\ ) and eq. ( p7l3| ) we conclude that 



Ul d ■ Uab ■ Uod — LU n da Uab 



U$a ■ Uab ■ U d0 = U^Uab 



(2.18) 
(2.19) 



i.e. U a Q (and their combinations) are invariant under 
cyclic permutations of states, the Uo a are invariant under 
cyclic permutations of phases. 

The determinant of U a b can only take on the values 
±1, 



detlWt-l)^"- 1 ) , 



(2.20) 



and the eigenvalues, consecutively numbered by k, all lie 
on the unit circle in the complex plane 



C. Complex coherence vector 

The general state of a quantum-mechanical system is 
specified by its density-operator p — ^ijPijPij with 

Pij = tr I pP}j I defining the respective density matrix. 

In many areas of physics, especially for finite-dimensional 
state spaces, the description of p using a coherence vec- 
tor (as the set of expectation values of the underlying 
SU(n)-operator basis) has shown great power because 
of its intuitive, almost "classical" behaviour. For two- 
level systems the coherence vector lives in an ordinary 
3-dimcnsional space. In general, the SU(ro) basis, which 
is hermitian, leads to a (n 2 — l)-dimcnsional vector con- 
sisting of real values, while the unitary operators studied 
here give complex elements. 

For the operators U a b we will henceforth use a single 
index notation Ui with i S [0,n 2 — 1], interpreting ab as 
the n-adic numb er rep resenting i. 

Following eq. ([T 121), we can expand the density oper- 
ator p as 



1 n -1 



UiUi 



(2.26) 



i=0 



and collect the coefficients to define the complex co- 
herence vector 



u := {ui , < i < n 2 - 1} 
u l :=(ul)=tv{lllp} EC 



(2.27) 
(2.28) 



where we deliberately excluded uq since it always 
equa ls 1 because of Uq — 1. From the symmetry (cf. 
eq. ( p3b ) 



Afe — LOn 



bl+ab- 



(-1) 



ab{n— 1) bl 



(2.21) 



From the cyclic symmetry the commutation properties 
are found to be 



U a b, U c d 



Uab,Ul 



(UJ U C ± W° ) Ug + c , b+d , 

, ,cdj , — bc _i_ —ad\fr 

V^n ± w n ) U a-c,b-d , 



and the structure constants f a b,cd,ef, defined by 



u ab ,u cd 



n-1 

^ > fab,cd,ef ^e/ j 
e,/=0 



follow as 

fab,cd,ef 



-be , ,—ad\ x x 



J e,a+cOf,b+d 



(2.22) 
(2.23) 

(2.24) 
(2.25) 



These structure constants are much simpler than those 
for SU(n)! Relations to other basis sets are summarized 
in Appendix B. 



Uab 



* ab 



(2.29) 



we conclude that there are n 2 — 1 independent real 
parameters forming the complex coherence vector. The 
length of u can be related to the trace of the density 
operator p 



|u| 2 = ntr{p 2 } - 1 , 



(2.30) 



thus giving a simple criterion to distinguish between 
pure (tr{p 2 } = 1) and mixed (tr{p 2 } < 1) states: 



< lul < n- 1 



pure state 
mixed state 



(2.31) 
(2.32) 



The convenience of the well known Bloch vector for- 
malism comes from the simple motion performed by the 
vector under unitary transformations: Since the length 



3 



|u| is preserved, any unitary time evolution operator U(i) 

just causes a rotation of the coherence vector. If the den- p = 
sity matrix evolves like p(t) = U {t) p(0)U {t)\ one can j-[ >, 

easily show that the motion of u is fj,=i 



N 



1 / N % " 



fj,<u i,j 



/n=l i=l 



u(t) = T(t) u(0) 
with rotation matrix T(t) M, 



(2.33) 



+E E-Sr'^ 




(12...N)f T (12...N) 
"i.7 A.../ u ijk...l 



T tJ (t) := ~tx{U]U(tyulu(t)} 



(2.34) 



Using the Liouvillc equation one can replace the time 
evolution operator by the Hamiltonian H{t) leading to a 
differential equation similar to the Bloch equations, 



u = Q u 



— tr{H(t) 
run 



(2.35) 



As expected, Q with its properties Qij = — f2*,j and 
tr{Oy} = describes a rotation in complex vector-space. 



III. CLUSTER-OPERATOR BASIS 

A. Definition and properties 

Up to now we have restricted ourselves to a single sys- 
tem, respective one node of a quantum network. To de- 
scribe a network of N subsystems, N > 1, looking at each 
node separately is not enough since correlations between 
nodes emerge - so the concept of clusters has to be intro- 
duced. A cluster operator acting on m particles is build 

out of m one particle operators tj^\ i ^ 0, where greek 
indices label different nodes of the network. All cluster 
operators 



:= 1 ( 



1 O Ui (g)l (g) . . . (g) 1 
node p 



node p 



node 



(3.1) 
i (3.2) 



etc. 



are unitary and orthonormal, and together span the 
complete Liouville-space of the quantum network. Act- 
ing on only m nodes out of iV means leaving the others 
unaffected by choosing Uo = 1 for them (Thus m is the 
number of indices unequal zero). The general de- 

composition of p (like of any other operator A) for a N 
particle quantum network with node v being a n„-level 
system reads 



(3.3) 



with the index-selective expectation values 



(pv...a) 

Hj...k 



= tr \p ■ U % 



ij...k 



} = (Ufk^) ■ (3-4) 



The local coherence vectors discussed in Sect. [IB 



just turn out to be the correlation tensors of first or- 
der, m = 1. If the network is in a product state, all 
correlation tensors of higher order factor into a product 
of local coherence vectors, 



(m) W 
uf'vq 



(3.5) 



i.e. these states are completely determined by local prop- 
erties. In case of a single subsystem, the length of its 
coherence vector ^\ \ui\ can be identified with a scalar 
constant of motion under unitary evaluation. For higher 
order correlation tensors, the concept of cluster sums 
leads to new invariants. For any m particle cluster, the 
cluster sum is defined as the weight of the tensor (cf . (|] ) 



y(liv...<j) ,_ (HV...CT) 




(3.6) 



with special cases Yq = 1 and = |u^| . For a 
network with N nodes, these cluster sums give 2 N scalar 
invariants under (products of) local unitary transforma- 
tions. All the cluster sums together are subject to the 
sum rule 



tr 



JV N 
fj,= l fi=l n<u 



(3.7) 



E y* 

H<v<o 



Y, 



(123. ..iV) 



N 



For a pure state, tr {p 2 } = 1 and for n M ~ n the left 
hand side is n N . According to eq. (3.5), any m-cluster 
sum of a product state factorizes into its 1-cluster-sum 
components. A given cluster of size m can then be tested 
as being in a non-product-state, if there is some partition 
into smaller clusters, the cluster-sum product of which is 



smaller than Y m , e.g. 



< Y 



(12) 



The cluster 



4 



sums can be related to the " purity" of a to particle clus- 
ter by interpreting the cluster as a n m -level system with 
coherence vector u m . (For simplicity we assume n M = n.) 
A normalized " purity factor" can now be defined from the 
respective coherence vector length as (m < N) 



'tr {Pciuster} 1 

n m - 1 



(3.8) 



i / m 

\H=1 li<v 



characterizi ng th e purity of a m-node cluster on the 
scale (cf. eq. fl2.32|)) 



<p„ 

Pm = 
Pm = 1 



< 1 



(3.9) 



maximal mixed TO-cluster 
pure TO-cluster. 



Alternatively, the purity could be characterized by the 
respective cluster-entropy [||; p m has the advantage of 
being a simple algebraic function of the expectation val- 
ues (matrix-elements of the reduced density matrix) . 



IV. COLLECTIVE OPERATOR-BASIS FOR 
TWO-LEVEL SUBSYSTEMS 

A. Definition and properties 

For now we will stick to quantum networks build out 
of N two level systems, knowing that a generalization to 
n-level systems is straight forward. First one needs to 
specify which single particle operator basis is used, e.g. 
uqi = a x , i>Qi = fry, w = a z (which is hermitian and 
unitary) or a± = a x ±ia y , a z or any other complete basis. 

We then relabel the cluster operators U^ 2 " k N ' 1 as C , Q( a 7iP , 
meaning an operator of dimension 2 2N , where a, (3 and 
7 specify the multiplicity of a x , a y , and a z , respectively 
(a + (3 + 7 < N). Index p specifies a permutation of 
these among the N subsystems. The number of such 
permutations and hence the index range for p £ [0, — 1] 



fi(a,)8,7) 



Nl 



a\f3\-y\{N - a - /3 - 7)! 



(4.1) 



These operators C a p-y. p , again, span the whole 
Liouville-space by defining all subsystem specific prop- 
erties. We now go on to define collective operators E 
by 



E, 



£7—1 

a/37,6 — '''' 
p=0 



pbs~i 
W f2 ( - / a/37,p ) 



(4.2) 



To ensure that all subsystems are treated on equal foot- 
ing, the sum extends over all permutations p, weighted 
only with pure phase fators, where b £ [0, fl — 1] labels 
the phase shift between "neighbouring" p. (Here the 
numbering of permutations is a matter of choice and the 
phase has no physical meaning.) The set of collective 
operators is orthonormal 



1 



tr {^a^.b-Ev^/y,;,/} — 5aa>5pl3<5 11 <<!)bb' (4.3) 



n2 N 

and complete, so that the density operator p of the 
network can be decomposed as 

P = E E E *fi 

{af3j} b 



y.{3'j 1 b-&a{3'y,b •> 



(4.4) 



E a f3y,b — 



(4.5) 



The expectation values E a /3j t b are collective in the 
sense that they do not refer to specific sub system -indices 
(as opposed to the £7-operators, cf. eq. (3.1,3.2)). The 
inverse transformation is given by 

£1-1 

V* / ,(p-po) b _ cin 

6=0 

(4.6) 



Q-l 

E W n C >E aP-l,b 
6=0 



n-i 

E 

p=0 



where we have made use of eq. (2.9). For example, the 
selective operator Cioo.o = <8>1<8>. • • could be written 



a-i 



E 



100,6- 



Note that the operators E are, in 



general, no longer unitary; for 6 = they are hermitian 
if the operators C a f3~ (tP are. 



B. Alternative sets 

The operators E still distinguish between all three ba- 
sic operators. Further reductions are possible: one such 
variant is (based on cr + , <r_, a z ) 



F zn , b -EE w n{z a y 



0(2,7) a-@=z 
P=0 a, 

which has been considered for nuclear spin- networks. In 
NMR the F z -terms in a Hamiltonian are said to induce 
a total of z quantum-" flip-flops" 0; an _F z -term entering 
the density-operator p describes the respective coherence 
order \z\. 

Another possibility is 



O(m) a+/3+f=m 

E E 



W n(m)^«/37,P ( 4 - 8 ) 



where p are permutations of the to operators (of any 
type) on the N subsystems, to = 0,1, ... ,N. Also this 
set is still complete. 
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V. APPLICATIONS 

A. Commuting operator sets and generalized cat 
basis 

1. Commuting sets of cluster operators 

Cluster sums and purity factors can be used to clas- 
sify states by their non-local properties. Highly entan- 
gled states tend to share correlations among all nodes 
rather than between only a few. Bell states |L(J {N = 2) 
are perfect in this sense, because they are in a totally 
mixed state (Yi = Pi = 0) locally and pure otherwise 
(Y 2 = 3,p 2 = !)• Generalizing this point, one may define 
"highly entangled states" as states that have a cluster 
sum distribution with a strong focus on multi-particle 
correlations. This definition does not give us a quantita- 
tive measure for multi-particle entanglement, but it can 
assist us in the search for new states with correlations of 
high order. Using the unitary operators defined above, 
we can give a constructive way of finding such states. We 
restrict ourselves to = n. 

Lemma 1 

Let \tp) be an eigenstate of n N completely commuting 
unitary operators out of which there are qjy 7V-cluster- 
opererators. Then its highest order cluster sum is given 
by Y N = q N . 

Note that the modulus of any eigenvalue of any cluster- 
operator U-j 1 "'"' is exactly f; the q N commuting TV- 
cluster-operators thus imply Yjv = qN provided all 
the non-commuting A-cluster-operators contribute zero. 
This must be the case, however, if there is a total of 



commuting operators exploiting the sum rule eq. (|3.7|). 

A set of completely commuting operators is a set in 
which each operator commutes with all others. This 
lemma reduces the problem of finding a state with max- 
imum cluster sum Yjy to the problem of finding a max- 
imum set of commuting A-cluster-operators. Before we 
can make statements about cluster operators, the com- 
muting relations of single particle unitary operators Ui 



need to be examined. Since 
0, one can show: 



ui b ,u cd 



= iff ad — be — 



Lemma 2. (based on theory of congruency classes) 
Each operator U a b commutes with ngcd(a, 6, n) other ba- 
sis operators and to each operator U a b there exists a triv- 
ial set of n completely commuting basis operators, (gcd 
= greatest common divisor) 

One of the n commuting operators in the completely 
commuting set is Uq = 1, so the maximal number of 
commuting one particle operators is n — 1. 

For a network of N particles with n-levels each, the 



properties of unitary cluster operators lead to the com- 
mutator relation 



U ai b! <S> U a2b2 <g> U Cldl ® U C2d2 <g> . 

a\di - bxci + a 2 d 2 - b 2 c 2 



= 
= 



(5.1) 



and from this we find the following statements: 



A There exists a set of (n — 1)^ completely commut- 
ing cluster operators of size N. This set can be 
constructed from all combinations of n — 1 com- 
pletely commuting one particle cluster operators. 



B There is a set of 

(n 2 

(n 2 - 1)— (n-1), 



1)- 



N even (5.2) 
N odd (5.3) 



completely commuting ./V particle cluster operators 
of the form 



Uaibi ® U biai <8> U a2 b 2 (g) U b2 a 2 -- 



(5.4) 



C There can be larger sets of co mple tely commuting 
operators than given by eq. (5.3). However, no 
general constructive method to find them is known 
to us so far. For n > 2 we restrict ourselves to 
N < 3. 

D The maximum order of any such set is constrained 
by n N - 1. This maximum value is typically not 
reached for N > 2. 

The results of A, B, C and D are summarized in Table 





n = 


= 2 












n - 


= 3 










N 


A 


B 


C 


D 


Cat 


N 


A 


B 


C 


D 


Cat 


1 


1 


1 


1 


1 


1 




1 


2 


2 


2 


2 


2 


2 


1 


3 


3 


3 


3 




2 


4 


8 


8 


8 


8 


3 


1 


3 


4 


7 


4 




3 


8 


16 


20 


26 


20 


4 


1 


9 


9 


15 


9 




4 


16 


64 


? 


80 


60 


5 


1 


9 


16 


31 


16 




5 


32 


128 


? 


242 


172 


6 


1 


27 


33 


63 


33 




6 


64 


512 


? 


728 


508 



n = 4 



A 



3 
9 

27 
81 
243 



B 



3 

15 

45 

175 

525 



C 



D 



3 

15 

63 

255 

1023 



Cat 



3 

15 

54 

213 

828 



TABLE I. Cluster sums Yn for eigenstates corresponding 
to sets of completely commuting N particle cluster operators 
calculated by methods A, B and C. Column D gives the upper 
limit and colum Cat states the cl uster s um Yjv for generalized 
cat states as introduced in Sect. VA2. 



G 



2. Generalized cat states 



As a first example let us look at the well known case 
of two spin 1/2 particles. The maximum order for com- 
pletely commuting sets of 2-cluster-operators is 3 and 
there are 6 such sets: 



(5.5) 
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The eigenstates to the first set form the Bell basis 
1 00) ± |11), |01) ± 1 10), the second set has i |00) ± |11), 
i |01)±|10) as eigenstates and the third gives |00) + |01) + 
|10)-|11), |00) + |01)-|10)+|11), |00)-|01)+|10)+|11), 
- 1 00) + |01) + 1 10) + 1 11). The remaining sets have eigen- 
states similar to the ones given by the third set, except 
for a phase. All these states are totally mixed (purity 
factor pi — 0) locally, and pure in total (p2 = 1). From 
the table for the cluster sums, one can see that in the 
case of spin 1/2 particles (n — 2) there are no states 
with more entanglement than the generalized cat states 
(see next chapter). For n = 3-level systems this is also 
the case for up to 3 particles, but in the case of 4 particles 
there exists a state with a higher cluster sum Y4 = 64. 
It turns out that this state is the dyadic product of two 
Bell states (each with Y<z = 8). This indicates, that a cat 
stat e may not necessarily have the largest Yjv-value (cf. 
eq. (p0D ). 

We generalize the Bell basis to systems with an ar- 
bitrary number of N particles where each particle is 
a n-level subsystem. Our starting point is the state 

n — l N 

\Cat) := — !^ J27=o ® N) w hich has the nice property, 

v k—l 

that any subsystem of m particles has an entanglement of 
S(pm) = log 2 n, when measured in terms of the local von 
Neumann entropy. (p m is the reduced density-operator of 
the respective m-cluster.) So any cluster, regardless of its 
size m, (to < N) has the same entropy and furthermore 
this entropy is the maximum amount of entanglement 
a n-level system can have. In this sense this state can 
be considered as a maximal entangled state. Since en- 
tanglement properties are invariant under local unitary 
transformations 0, one can generate a complete basis 
set from |Cat) by applying the discrete set of unitary 
basis operators Ui. The result, 



3=0 



(5.6) 



is an explicit definition of the cat state basis. The index 
c = {c,i}, Ci € [0, n — 1] labels the n N states and u> n is 
given by eq. ( |2.8| ). Orthonormality and completeness 
read 



d (Cat \Cat) c - 5 cA , ^ \Cat) c c (Cat \ = 1. (5.7) 



In the low dimensional case of N = n = 2 the definition 
( |5.6| ) reduces to the bell basis and for N = 3, n = 2 the 
GHZ state i(|000) + 1 1 1 1) ) is one member of such a set. 
Further examples for cat-states are given in Appendix C. 

The characterization of these generalized cat states in 
terms of cluster sums shows (cf. Table B) 



F x = 0, 

yw= (n-l)" + (-!)"(»-!), 1 < m<N) 



(5.8) 
(5.9) 



y w = (n-l)n JV - i + 



N-l , (n-l) N + (-l) N (n-l) 



(5.10) 



while the puri ty fa ctor distribution for any such a cat 
state is (cf. eq. (3.8)) 



Pi = 0, 



m— 1 



Prn 
PN 



n" 
1 , 



(totaly mixed) (5-11) 
1 < to < N (5.12) 
(pure) . (5.13) 



The purity factor distribution for cat-states is shown 
in Fig. @. 

For n> 1 any cluster, regardless of its size to, (to < 
N) approaches the same purity factor p m — —. This 
behaviour is reminescent of that for the cluster-entropy. 

For N — > the cluster-sum Yjy almost exploits the 
sum rule (3.7): 



n N 



1 



(5.14) 



Cat states thus become even more "non-classical" with 
increasing TV (and n). 

B. Collective invariants and eigenstates 

We have already briefly discussed invariants under gen- 
eral unitary transformations (tr{/5 2 }) and under local 
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unitary transformations (the cluster sums). Let our 
quantum network be described by the Hamiltonian H, 
with the spectral representation 



H — ^ Ekf>k 

k 



(5.15) 

where pk = \4>k) (<^fe|- Then the expectation- values 



C J =tr{ppj} = (^\p\h) (5.16) 

constrained by Y2j Cj — 1 are additional invariants 
of motion under the unitary transformation generated 

by H. This follows from the fact, that i?i^tr|pA| = 
tr I [H, p] _ A I = tr -I [A, H] _ p j which vanishes if A com- 
mutes with H. Note that \(f>j) has, in general, not product 
form. We consider 3 examples: 

i. N = n = 2 network with Forster-interaction Cp 
The Hamiltonian in terms of the collective opera- 



tors (cf. Sect. IV A) 



H = —r- £001,0 H ^ — (£200,0 + £020,0) (5-17) 



has the eigenstates 

\<M = loo) 



\<fi±) =^(|01)±|10» 

I01> = 111) 



which can be rewritten as 
1 



Pa 



(i + £, 



002,0 



£001. c 



P± — j(l — £q02,0 ± (£200,0 + ^020,0)) 



Pi 



-(! + £, 



002,0 



E, 



001.0) 



(5.18) 
(5.19) 
(5.20) 

(5.21) 
(5.22) 
(5.23) 



These are all permutation-symmetric density- 
operators; their expectation values Cj according to 
eq. (5.16) lead to three independent collective in- 



variants 



£001,0 


= const. 


(5.24) 


£-002,0 


= const. 


(5.25) 


£200,0 + £-020,0 


= const. 


(5.26) 



ii. N — n = 2 with renormalization interaction Cj^ 
The Hamiltonian 



H= - (J 1 ' W 2) )£ooi,o 



(5.27) 



-i- h (, ,(!) ,,( 2 )^p 1 hC R ft 

+ — & J-C'OOl,! H 2 — fc 002,0 



has the product states as eigenstates. The collec- 
tive invariants are in this case 



£001,0 


= const. 


£-002,0 


= const. 


£■001,1 


= const. 



(5.28) 
(5.29) 
(5.30) 



iii. N = n = 2 with collective stimulation 
The Hamiltonian given by 



H — - g £100,0 



— 5 £001,0 



(5.31) 



describes a non-interacting spin-pair driven by a co- 
herent electromagnetic field (strength g, detuning 
5) in rotating wave approximation. Its invariants 
turn out to be 



^£001,0 + gE\ 



const. 



(5.32) 



S 2 £002,0 + 2g5E 101 . + g 2 E 20 o,o = const. (5.33) 

4(?V + <5 2 )£ooi,i - 5 4 £oo2,o 
-9 2 {9 2 + S 2 )E 020 ,o + 4g6(g 2 + S 2 )E W0A 

+2g 3 SE wlfi - g 2 S 2 E 2 oo,o = const. (5.34) 



Contrary to example ii, the _ff-model ( 5.31 ) is 
permutation-symmetric, while the invariant (5.34) con- 
tains expectation values with b =/= 0. The latter can be 
unequal zero only if at some previous stage of the prepa- 
ration the permutation symmetry had been broken (cf. 
Sect. p|. 



C. Permutation symmetry 

A complete description of large quantum networks 
turns out to be virtually impossible, as the Hilbert space 
dimension grows exponentially with the number of parti- 
cles used. However, symmetry can reduce the complexity 
of the total system by a significant amount: For a per- 
fect permutation-symmetric system we will show that the 
reduction in the number of parameters needed is so enor- 
mous, that a polynomial increasing number of parame- 
ters is enough to describe networks of arbitrary size. 

Within such a network no subsystem can be distin- 
guished, neither in preparation nor in detection. Two 
possible setups can be thought of: Fundamental indistin- 
guishable and operational indistinguishable subsystems. 
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A typical setup of the former type would be several elec- 
trons in a box. As the electrons are fermions and their 
location is not fixed, there is no way to act independently 
on a specific electron. All controlling and measurement 
procedures act on the system as a whole. Therefore no 
information loss can occur when reducing the system- 
description to permutation-symmetric operators. 

Opposed to that, in the case of operational indistin- 
guishability, the design of the experimental setup is the 
source of reduction. A linear ion trap with a laser beam 
acting on the ions could be used as an example. Only if 
the beam waist is less than the spatial separation of the 
ions, the particles become distinguishable, so by control- 
ling the laser beam the experimentalist can choose which 
case he is in. 

Collective operators E a /3-y,b with 6 = (cf. Sect. IV A) 
are invariant under any arbitrary subsystem permutation 
P 



(5.35) 



For indistinguishable subsystems, these are the only op- 
erators allowed: Any control or measurement operator 
has to be part of the Liouville-space spanned by Eap^fi- 
The number of such operators is (for n = 2) 



(JV + l)(JV + 2)(jV + 3) 
6 



(5.36) 



and thus scales polynomially of order 0(N 3 ) with the 
number of subsystems involved. Therefore, such a highly 
symmetric system constitutes a class of reduced complex- 
ity as compared with a general quantum network of the 
same size N. This is somewhat surprising as indistin- 
guishability might have been expected to enhance non- 
classical features. 

Alternatively, if all expectation-values happen to be 
permutation-symmetric (no "structure"), then 



(Ca/?7,p) — (Ca/37) 



for all jp-perm utat ions, and by recalling eq. 
definition eq. (4.2) leads to 



(E a p ly b) — Q(C a fjy}5ob 



(5.37) 
|) the 

(5.38) 



i.e. we need to consider only the £ a /3 7 ,o operators. This 
would hold, correspondingly, for the operators F Zil ^ ac- 



cording to eq. (4.7) (or G m ,b according to eq. (4.8)) pro- 
vided only z, 7 (or m) could indeed be distinguished. The 
number of remaining operators is then further reduced: 



For F z „ fi this number is E 7 =o( 2 ( iV + = + 
for G m> o it is simply N + 1. 



1. Symmetry classes 

Permutation symmetry defines a kind of operational 
indistinguishability between the subsystems of a given 
network. This symmetry alone, however, would allow 
for more symmetry classes than realized in nature by 
the fundamentally indistinguishable particles, Fermions 
and Bosons, respectively. The spin-statistics- relation go- 
ing back to Pauli jy],]^ m ight be relaxed, if the parti- 
cles or subsystems are localized in different areas of real 
space |fl2|j . The location index would, in principle, ren- 
der these subsystems distinguishable; however, for the 
following we assume that the actual operators describ- 
ing the network and its coupling to the outside world 
are still permutation-symmetric, so that the correspond- 
ing super-selection rules apply, as will be discussed be- 
low. This would imply, e.g. that electrons, localized in 
different semiconductor quantum dots, could live in the 
state-subspace of Bose-symmetry (or of any other "para- 
boson" |13| symmetry class as well; cf. Table ID), pro- 
vided we are able to prepare such states from some stan- 
dard initial state: Directed transi ent symmetry breaking 
could do this job (cf. Sect. |VE4| ). 

Any symmetry type can be characterized by a Young 
diagram and is equivalent to a irreducible representation 
of the permutation group Sn spanned by the basis vec- 



tors. As an example, Table III shows all possible sym- 
metries for a N = 4 particle system. Since a definite an- 
gular momentum quantum number j is assigned to every 
Young diagram, there are 2j + 1 states of equal symmetry 
type but different energy (i.e. different configuration). If 
the system is subjected to permutation-symmetric oper- 
ators only, the super-selection rules prohibit any transi- 
tion between different Young tables. In case the system 
under consideration is in a specific state with angular 
momentum j, the state space reached by applying any 
collective operator Eijk.o is of dimension 2j + 1. There- 
fore (2j + l) 2 parameters (expectation values) are needed 
to describe the system. The total number of parameters, 

E£o( 2 J + !) 2 = U N + X )( N + 2 )( N + 3 ) is exactly the 
number £0 of collective operators Eijkfi (see eq. ( 5.36 )). 



2. Structure and Hamilton-models 

For fundamentally indistinguishable particles, the sub- 
system index fj, has no physical meaning. However, it 
may happen that a specific property is not only a good 
quantum number but a unique constant of motion for any 
subsystem. Localization in real space is a pertinent ex- 
ample; the subsystem-index [i is then mapped onto a spa- 
tial position-index R^. By this position any subsystem 
becomes distinguishable, in principle. The phases (en- 
tering the collective operators E) get a physical meaning 
("wave length"). An "operational" indistinguishability 
remains, if the Hamiltonian describing the network does 
still contain collective-operators of the type -B Q( 3 7 ,o only. 







Typical Hamilton- models include m — 1,2 -particle 
operators. The structure tends to break permutation- 
symmetry for the localized states as the coupling usually 
depends on the distance \R^ — R v \. This partial selec- 
tivity can be described as a perturbation via collective 
operators E a fj 1 ^ with b ^ 0. 

If all pertinent distances could be made equal, the 
breaking of the permutation symmetry would go to zero. 
For N > D + 1 (in D dimensions) however, the interac- 
tion distances cannot all be the same. A partial remedy 
is the introduction of a " quantum bus" : In this case the 
nodes do not interact directly but only indirectly via a 
common degree of freedom. This degree of freedom could 
be a central spin, but typically is implemented as a col- 
lective mode (like the phonon mode of a cold ion trap 



D. Symmetry breaking and irreversibility 

Constrained operations will lead, quantum mechani- 
cally, to selection rules accompanied by a tremendous 
reduction of the state space available to the system dy- 
namics starting from a given initial state. This situation 
needs to be distinguished from lack of control (measure- 
ment data) implying lack of information (entropy S > 0). 
In so far as this lack of control refers to expectation val- 
ues 6^0 (which would be absent for strict permutation 
symmetry), "uncontrolled" symmetry breaking may be 
said to lead to an ensemble description and irreversibil- 
ity: The true state space has been reduced to the smaller 
one defined by the assumed permutation symmetry. In 
this case, however, the " decoherence" does not reflect the 
influence of an external bath but is rather of "internal" 
origin. Phenomenologically ony may try to model these 
effects, as usual, via some decoherence times; however, 
it is not clear yet under what conditions such a proce- 
dure would be appropriate. Clearly it should, at most, 
work for sufficiently large networks and only to the ex- 
tent that the reduction is really substantial (cf. ref . [jl5| ) . 
The decoherence time would then, in turn, allow to assess 
symmetry breaking effects in a global way. 

In the following, we will restrict ourselves to unitary 
transformations. 



E. Special unitary transformations 

1. Cyclic permutations: generalized echoes 
The pulse-like manipulation of Hamiltonians, 



H(t) = H 3 , for tj<t< tj+i , j = 0, 1, 2, . 



(5.39) 



U 3 {At 3 



At, 



t 



3+1 



u 



(5.40) 



at will, has become a popular approach to quantum 
computation || . One basic operation is the "halting op- 
eration", Uj(Atj) = 1, which would require Hj = 0. 
This condition cannot be realized, in general. The well- 
known spin-echoes [fl7f are formally based on time rever- 
sal, thus undoing a unitary evolution over a time period 
At/2 within the next period At/ 2: 



U 3 {-At/2) = (j]{At/2) 



(5.41) 



This time-reversal cannot be implemented either; how- 
ever, its effect can be simulated in any discrete Hilbert- 
space of dimension n by means of the cyclic perm utation 
operations U n -i Oj &s introduced in section II B. 



Let us consider the time-independent Hamiltonian H . 
It is convenient to require tr = so that 



(5.42) 



with Ek denoting the eigenvalues of H. Then the effect 
of the unitary time evolution generated by H 



|V(Ai)) = U 6 (At) 1-0(0)) (5.43) 
can be suppressed for any given At by the iteration 



U eS (At) := [U n - li0 U 6 (At/n) 



(5.44) 



Here we have assumed f or s implicity that the appli- 
cation of U n -ifl (cf. eq. (2.8)) does not consume any 
additional time. U n -i^ generates cyclic permutations 
between all the states Ek of the spectrum. Any ini- 
tial eigenstate thus visits all the other eigenstates f or the 
same time At/n so that, due to the constraint eq. ( 5.42 ) 
the total acquired phase adds up to zero. This invariancc 
property holds for any initial state \ip(0)}, as such a state 
can be written as a superposition of eigenstates, and any 
eigenstate returns to it s ini tial phase separately. 

This mapping, eq. ( |5.44| ), can be repeated to get the 
stroboscopic invariance (cf. ref. p6[) 



U eS (mAt) = (u n - 1>0 U (At/n 



= 1 



(5.45) 



to shape unitary evolutions 



where any initial state is periodically recovered. The 
cyclic permutation U n —i t o can be decomposed into a se- 
quence of n — 1 state-selective 7r-pulses; for n = 3, e.g.: 
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#2,0 = Pox + P12 + P20 (5.46) 
= (-P00 + P12 + P21) ■ {P01 + P10 + P22) 

The total number of such 7r-pulses for U e ff(mAt) is 
n(n— l)m. The conventional spin-echo obtains for n = 2, 
m = 1, and an initial state | -0 (0)} generated from the 
ground state by means of a 7r/2-pulse. The second U\fi- 
flip is then usually omitted. 

2. Selective control 



N\s 





1 


2 


3 


4 


5 


6 


7 8 . . . 


1 





1 














2 


00 


01 


11 


10 










3 


000 


001 


011 


010 


110 


111 


101 


100 


4 


0000 


0001 


0011 


0010 


0110 


0111 


0101 


0100 1100 ... 



TABLE II. Circular sequences of product states with Ham- 
ming distance 1 for various network sizes N. These sequences 
can be constructed recursively: Sequence N + 1 is obtained 
from that of N by repeating each member and supplementing 
it from the right by 0, 1 (1, 0), respectively. (There are other 
sequences with the same property.) 



In principle, the results obtained in Sect. VE1 may 



also be used for a network of N subsystems of dimen- 
sion n each. The condition is, that one is able to 
implement the respective cyclic operator: For a net- 
work with (distinguishable) non-interacting subsystems 
we simply have to replace U n -i by the product-operator 
r(i) 



Q — 



U n-1,0- 



In the case of identical subsystems, interactions may 
lift the degeneracy of the network-eigenstates. We con- 
sider the example of a network of 2-level-subsystems and 
an interaction of the type 



H' 



h c »' x 



R A 3 



(5.47) 



This interaction modifies the eigen-spectrum (allowing 
for selectivity) while leaving the eigenstates unchanged 
(i.e. as p roduct-states), \p(N) . . . p(l)); p(p) =0,1 (cf. 
eq. ( 5.27 )). Single-particle transitions are allowed only 
between states which differ at one position ("Hamming 
distance" 1). A convenient numbering, s, of the 2 N states 
allowing for cyclic permutation is given in Table |l[ For a 
large network N, the echo becomes difficult to implement 
practically, as it requires ~ 2 2N short but at the same 
time frequency-selective 7r-pulses per time period At. In 
general, interactions will also modify the eigenstates (i.e. 
lead to a non-product for m). A si mple example for N = 
n = 2 is provided by eq. (5.18) - (5.2C). 

It has been proposed ]18| to exploit some approxi- 
mate stroboscopic invariance also for open systems un- 
der the condition At <C r c , where r c is a typical time- 
scale for dissipation. However, such short pulses will 
easily come in conflict with their required selectivity in 
frequency-space and thus tend to restrict practical im- 
plementations of their proposal to small networks only. 
Furthermore, it is unlikely that the additional control 
interactions could work without inducing uncontrollable 
features on their own. Supression of decoherence would 
amount to increase the pertinent decoherence-time and 
-length to macroscopic dimensions (which typically hap- 
pens for very special and restricted state-spaces only, like 
in superconductivity or Bose-Einstein condensation). 



3. Collective control 

We consider the unitary evolution induced by a "col- 
lective" and permutation-symmmetric Hamiltonian like 

H = haE m00fi = ha ^ ■ --^ (5.48) 

This H generates the unitary time evolution 



U(t) 



-iatE„ 



(5.49) 

= J] (cos(ot)i-iBin(at)^M...d-«)). 
n<v...<(; 



In particular for m = 1 and at — -| on finds U* — 
(-i) Enoo,o, while for at = f , U* cx Y^=o (-*)"^oo,o- 
For m = 2 and at = \ , we have Us. = (— i)~£jvoo,o 

JV — 1 A 

for N even and U * = i 1 otherwise; for N even and 
time at — j, eq. ( 5.49 ) simplifies to 



N 

Utl all i-EWochO , for — even (5.50) 

N 

Us. oc 1 - iE N00 fi , for — odd (5.51) 

Applied to the ground state |00...0), thi s Us thus 
creates the Af-particle cat state (cf. Sect. VA2 and 
Appendix C) 



Us |00...0) oc |00...0) ±i|ll...l) 



(5.52) 



in a single step! Specific examples for such collective 
control scenarios have been discussed in Refs. 1191- 



4- Interplay between selective and collective interactions 

Selective interactions violate the selection rules im- 
plicit in permutation-symmetric interactions. This qual- 
itatively different dynamical behaviour can be exploited 
to implement specific functionalities: One possibility is 
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to address different symmetry-ciasses via seiective cou- 
pling (i.e. controlled symmetry breaking); another pos- 
sibility is to suppress transition due to permutation- 
symmetric interactions by using states of different sym- 
metry, which, nevertheless, could all be prepared selec- 
tively ("decoherence-free subspace"). 

To be specific, let us consider a network of N pseudo- 
spins without mutual interactions, but in the presence of 
a quantum bus (i.e. a collective mode to which all spins 
are coupled in the same way). We further assume that 
the coupling to a larger field can be made at will either 
selective or collective. We start (for N — 4, n = 2) with 
the permutation-symmetric ground state |0000), (j = 2). 
By applying selective laser pulses and exploiting the cou- 
pling to the quantum-bus we can generate the EPR-state 
-7j(|01) - 1 10)) within the pair v = 1,2. The total 4- 
particle state then becomes a member of the symmetry 
class j = 1, see Table III. If permutation symmetry 



is restored now, transitions are possible only within the 
(2j + l)-dimensional subspace of this symmetry class (on 
a time-scale less than the decoherence time). This would 
allow to study finite systems of (operationally) indistin- 
guishable subsystems of a symmetry not realized in na- 
ture by fundamentally indistinguishable particles! It is 
straightforward to extend this scheme to N > 4. 

Alternatively, the selection rules can readily be ex- 
ploited as a means for stabilization: For this purpose 
we assume that the "unwanted" dynamics (coupling to 
the bath) is permutation-symmetric, while the control 
dynamics to be used is selective (i.e. not subject to the 
selection rules). For N = 4 we could take the three low- 
est energy states j — 2,1,0. Unitary dynamics within 
this subspace would (in the ideal case) not be perturbed 
by dissipation, if the bath was kept at zero temperature. 
(There are higher energy levels with the same symmetry, 
though.) 

It is thus preferable to use the multiplicity of the states 
j = (for N even) which would decouple from the bath 
exactly, as each of those states is the only member of 
its symmetry class (see Table |l|): There is nothing to 
connect to under the action of a permutation-symmetric 
coupling (cf. rcf p5|). Such schemes have been inves- 
tigated by a number of authors |2^-|30||. The resulting 
stabilization is limited by the fact that the symmetry 
selection rules will, in practice, not hold strictly. 



operators have been derived. The latter are adapted to 
situations in which individual subsystems cannot be se- 
lected. Classically this would imply a loss of informa- 
tion about the resulting state. Quantum-mechanically, 
this lack of control would rather give rise to entangle- 
ment - with super-selection rules tremendously reduc- 
ing the state-space available: Actually its dimension be- 
comes polynomial in N. Such states have been proposed 
as subspaces of reduced decoherence (with respect to 
permutation-symmetric bath interations). We have sug- 
gested that localized, artificial quantum networks, sub- 
jected to certain symmetry-breaking initialization steps, 
might live in state-symmetry-classes that otherwise can- 
not occur for fundamentally indistinguishable particles. 
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VI. SUMMARY AND CONCLUSIONS 



In this paper we have addressed composite systems 
consisting of N subsystems [i of respective dimension n^. 
We have been concerned with applications of (selective) 
cluster-operators and collective operators, both based on 
unitary local operator sets, and both being orthogonal 
and complete. The former are very useful in characteriz- 
ing entanglement of pure network states, in particular of 
cat-states, for any n and N. Commuting sets of cluster- 
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APPENDIX A: UNITARY OPERATOR BASIS FOR N = 3 
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APPENDIX B: RELATIONS BETWEEN OPERATOR SETS 

The relations between this unitary operator set |t/ a &j, the SU(n) generators {ujk, Vjk, w{\ and the projection 
operators j Pij j are given by: 



P ^ = \^n bk U^ ( B1 ) 



6=0 



U ab = ]T ^ b Pa+k, k (B2) 



fc=0 



C>ah= Yl \( UJ n 5 ^k+^8k ± a, 3 )u jk + (B3) 



2 

0<i<fc<n 
0<j<k<n 



-, 2 + 1 



0</<t-i V2(i + !)(/ + 2) 

1 n— 1 

fi i* = - E « 6fe 2z_ + w » M ) (B4) 



a,h=0 
. n-1 

%fe = - H ^« b?£ ^.J-fc _ < bl $a,k-j) Uab 
a,fc=0 

/ " n-1 ;+i 

V 6=0 9=1 
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APPENDIX C: EXAMPLES FOR GENERALIZED 
CAT STATES 

We give some (low dimensional) examples for general- 
ized cat states consisting of N subsystems with each one 
being a n level system. Together they form a orthonor- 
mal and complete basis of the underlying n N dimensional 
Hilbert space. 

N = 2: 



n = 2: 



\Cat) m 
\Cat) m 
\Cat) w 
\Cat) n 



^ (100) + |11)) 
^(|01> + |10» 
^(100) -|11)) 
^(|01>-|10» 



n = 3: 



\Cat) 00 = ^(|00) + |11) + |22)) 
|Cat> 01 = ^(|01> + |12> + |20>) 
|Co*>02 = ^(|02> + |10) + |21)) 

\Cat) w = ^ (|00) + ei™ |11) + e-i" |22) 

\Cat) u = ^ (J01) + ei™ |12) + e"§™ |20) 

|Cot) 12 - ^ (|02) + e§™ |10) + e-i" |21) 

|Cai) 20 



| Cat), 



I ^u, 1-/22 

N = 3: 

n = 2: 



= ^(|00) + e- 


"i™ |11) + ei™ 


= 7l(l 01 >- 


fe" 


"§ 7ri |12) + 65^ 


= ^(|02) + e- 


"i™ |10) + ei™ 


K? a *)ooo = 


i 


(|000) + |111)) 


C a *)ooi = 


i 


(|001) + |110)) 


K? a *)oio = 


1 

V2 


(1010) + |101» 


|Cat) 011 = 


1 
72 


(|011> + |100» 


\Cat) W0 = 


1 


(|000> - |111» 


\Cat) wl = 


1 

71 


(1001) -|110)) 



|20) 



|Cot) 110 = ^(|010)-|101)) 
|Cat) m = ^=(|011)-|100)) 



The cat-states N = n = 2 can entirely be expressed in 
terms of 2-particle collective operators: 



Poo = 


-(1 + £200,0 


- £020,0 


+ £002,0) 


(Cl) 


Poi = 


-(i + £200,0 


+ £020,0 


- £002,0) 


(C2) 


pw = 


-(i - £200,0 


+ £020,0 


+ £002,0) 


(C3) 


hi = 


-(i — £200,0 


- £020,0 


— £002,0) 


(C4) 



These are all permutation-symmetric and involve N = 
2-particle operators only. For N > n = 2 collec- 
tive operators appear also of order < N: For N = 3, 
|Cat) 000 = ^(|000) + |111», e.g., one finds 



Pooo — 23 (1 + -^300,0 — £120,0 + £002,0) (C5) 

The last term is responsible for the surviving 2-particle 
properties. In general, the cat-states are no longer 
permutation-symmetric: A pertinent example is 



1 ~ 1 - 2 - 

POOl =2~3 (1 — g£002,0 + ^£002,1 



g-£l20,0 



3 



120,1 



2 - 

g£002,2 (C6) 

2 _-7T - 

— e 1 3 £120,2) • 
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I III 

FIG. 1. m-cluster purity factor p m , as defined in eq. (|3.S|), for cat-states depending on the subsystem dimension n and the 
cluster size m for m < TV. The total cat-state is a pure state, so pn always equals 1. 
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TABLE III. Symmetry classes for N = 4 subsystems. The configuration reflects the number of spins up or down and 
therefore the energy level m. The angular momentum quantum number j defines the symmetry class with respect to permu- 
tation. The basis vectors are uniquely labeled by their Young tableau together with the magnetic quantum number m. For 
permutation-symmetric operators, super-selection rules only allow transitions between states with the same Young tableau. 
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